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ABSTRACT 



Context. The diffuse interstellar medium (ISM) is an integral part of the evolution of the entire Galaxy. Metals are produced by stars 
and their abundances are the direct testimony of the history of stellar evolution. However, the interstellar dust composition is not well 
known and the total abundances are yet to be accurately determined. 

Aims. We probe ISM dust composition, total abundances, and abundance gradients through the study of interstellar absorption features 
in the high-resolution X-ray spectra of Galactic low-mass X-ray binaries (LMXBs). 

Methods. We use high-quality grating spectra of nine LMXBs taken with XMM-Newton. We measure the column densities of O, Ne, 
Mg, and Fe with an empirical model and estimate the Galactic abundance gradients. 

Results. The column densities of the neutral gas species are in agreement with those found in the literature. Solids are a significant 
reservoir of metals like oxygen and iron. Respectively, 15-25% and 65-90% of the total amount of Oi and Fei is found in dust. 
The dust amount and mixture seem to be consistent along all the lines-of-sight (LOS). Our estimates of abundance gradients and 
predictions of local interstellar abundances are in agreement with those measured at longer wavelengths. 

Conclusions. Our work shows that X-ray spectroscopy is a very powerful method to probe the ISM. For instance, on a large scale 
the ISM appears to be chemically homogeneous showing similar gas ionization ratios and dust mixtures. The agreement between the 
abundances of the ISM and the stellar objects suggests that the local Galaxy is also chemically homogeneous. 

Key words. ISM: abundances - ISM: dust, extinction - ISM: molecules - ISM: structure - X-rays: ISM 



1. Introduction 

The interstellar medium (ISM) is one of the most important com- 
ponent of galaxies because it influences their evolution: the ISM 
is enriched with heavy elements by stellar evolution, and it pro- 
vides the source of material for following generation of stars. 
There are several physical processes that alter the metallicity of 
the ISM. Stellar winds and supernovae expel part of the interstel- 
lar gas out of the Galactic disk, but gravity generally forces the 
gas to fall back through the process known as "Galactic foun- 
tain" (Shapiro & Field 1976). Gas accreted from smaller galax- 
ies, like the Magellanic Clouds, and the intergalactic medium in- 
creases the reservoir of low metallicity gas. The ISM also shows 
a complex structure consisting of phases at d ifferent equilibrium 
temperatures (for a review, see lDrainell201ll) . The cold phase is 
a blend of dust, molecules and almost neutral gas below 10 4 K. 
The warm and hot phases are mostly gaseous. The warm gas is 
weakly ionized, with a temperature of ~ 10 4 K. The hot gas is 
highly ionized, with temperatures of about 10 6 K. These three 
main phases are not entirely separated from each other. For in- 
stance, it is thought that a conductive cooling layer, revealed by 
Ovi, lies bet ween the warm and hot ionized phases (see e.g. 
lRichterll2006h . 

In the spectra of background sources the ISM produces red- 
dening and absorption lines. High resolution X-ray spectroscopy 
has become a powerful diagnostic tool for constraining the 
chemical and physical properties of the ISM. The K-shell transi- 
tions of carbon, nitrogen, oxygen, neon, and magnesium, and the 
L-shell transitions of iron fall in the soft X-ray energy band. The 



presence of different ionization states constrains the multi-phase 
temperature structure of the ISM. A summary of the results ob- 
tained b y X-ray spe c trosco py on the ISM in the last decades is 
given bv lPinto et alj d2010h. 

Briefly, Schattenburg & Canizaresl ([1986) first constrained 
interstellar oxygen in the X-ray band using the Einstein 
Observatory, but a thorough study of the ISM was possible 
only after the launch of the XMM-Newton and Chandra satel- 
lites, provided with high-resolution gratings. Absorption struc- 
ture due to ionized gas and dust was found around the inter- 
stellar oxyg en K-shell absorption edge in the spectra of sev- 
eral sources (iPaerels et al.ll200 It Ide" Vries et al.ll2003HJuett et all 
| 2004tlCostantini et alJl2005l .1201 2tlde Vries , & Costantinil l2009: 
iKaastra et aJ.ll2009fa iPinto et alj|2010ll2012allbl) . 

However, in the last decades, the most important discover- 
ies on the physics of the ISM phases have been obtained with 
UV, IR, optical, and radio spectroscopy. Dust is usually studied 
in IR, molecules in radio and IR, while the neutral and warm 
phases of the i nterstellar gas are also probed with UV data (see 
iFerrierelbOOll and references therein). The hot gas is the only 
ISM phase which is more often studied in X-rays. However, de- 
spite the very high resolution and signal-to-noise ratio at long 
wavelengths, both the chemistry and physics of the interstellar 
medium are still under debate. Elemental abundances, dust de- 
pletion factors and dust chemistry are not yet well determined 
and thus provide an open and interesting research field. The 
heavy elements such as oxygen (the most abundant one) and iron 
are produced in high-mass stars. Molecules and dust compounds 
are mainly formed in AGB stars and then grow in the diffuse 
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ISM (see e.g. Mat tsson & An dersen 2012j). Heavy elements di- 
rectly witness the history of the past stellar evolution. In the last 
years there have been many attempts to measure the interstellar 
abundances, but in many cases only lower limits were obtained 
due to the depletion of the elements from the gaseous phases into 
dust grains and to strong line saturation at long wavelengths. The 
dust depletion factors and especially the compounds mixture are 
not yet well known and thus the total ISM abundances are af- 
fected by strong systematic uncertainties. For instance, a signif- 
icant fraction of the oxygen which is bound in the solid phase 
is still unknown: silicates and carbonaceous compounds are not 
abundant enough to cover the amount of depleted oxygen. Ices 
may solve this problem, but they are hard to detect ( for a review 
about depletion factors in the ISM see Jenkins 2009). Therefore, 
our work focuses on an alternative method to determine ISM 
column densities, dust depletion factors, molecular compounds, 
total abundances (gas + dust), and abundance gradients through 
X-ray spectroscopy. 

In this work we report the detection and modeling of in- 
terstellar absorption lines and edges in the high-quality spec- 
tra of nine LMXBs: GS 1826-238, GX 9+9, GX 339-4, 
Aql X-l, SAX J1808.4-3658, Ser X-l, 4U 1254-690, 4U 
1636-536, and 4U 1735-444. The observations are taken 
with the XMM-Newton Reflection Grating Spectrometer (RGS, 
Iden Herder et al.l [20011) . see Table [1] for details. These sources 
are well suited for the analysis of the ISM because of their 
brightness and column densities A^h ~ 1 - 5 x 10 25 trT 2 (see 
Table [3), which are high enough to produce prominent O, Fe, 
and Ne absorption edges in the soft X-ray energy band. This pa- 
per provides a significant extension of the previous work on the 
ISM in the line-o f-sight (LOS) towards GS 1826-238 done by 
iPinto et al] d2010l) . We have used updated atomic and molecular 
database and extended the analysis to several LOS (see Fig.[TJ. 

Our analysis focuses on the 7 - 38 A first order spectra of 
the RGS detectors. We have not used the EPIC spectra of the 
observations because they have much lower spectral resolution 
and we were not interested in determining the source contin- 
uum at energies higher than the RGS band. Most of interstellar 
X-ray features are indeed found in the soft X-ray energy band. 
We performed the sp ectral analysis with SPEXQ version 2.03.03 
dKaastra et alj! 996). We s caled the elementa l abund ances to the 
proto-Solar abundances of lLodders & Palmel (120091) . We use the 
^-statistic throughout the paper and adopt lcr errors. 

The paper is organized as follows. In Sect. |2]we present the 
data. In Sect. [3] we report the spectral features that we analyze. 
In Sect. [4] we describe the models we use and the results of our 
analysis. The discussion and the comparison with previous work 
are given in Sect. [5] Conclusions are reported in Sect. |6] 

2. The data 

For this work we have used 22 archival spectra of nine sources 
taken with XMM-Newton (see Fig.[T]and Table[TJ. These spectra 
were taken in periods of high-flux state of the sources and each 
single spectrum shows high statistics. We have excluded in our 
analysis additional observations taken during low-flux states. For 
some sources we have just one good spectrum, while for other 
ones we obtained up to five high-quality spectra. The spectra 
taken on the same source have been fitted together in order to 
increase the statistics, but we have used a different approach be- 
tween stable and variable sources. When the persistent emission 
of a certain source was steady and the spectra were perfectly 

1 www.sron.nl/spex 



superimposable, we hav e stacked the spect ra according to the 
procedure described by Kaa stra et al.l (1201 ll) . Briefly, we have 
reduced the spectra with the Science Analysis System (SAS) ver- 
sion 11.0, obtained fluxed spectra for RGS 1 and 2 and stacked 
them with the SPEX task RGS- Huxcombine. We cannot stack 
the spectra of variable sources because this may introduce spuri- 
ous features near the interstellar edges. Thus, we simultaneously 
fit their spectra coupling the parameters of the interstellar ab- 
sorbers, but adopting a different continuum for each observation. 
This can be done with the SPEX task sectors. 

The LMXBs GS 1826-238, GX 9+9, Ser X-l, 
4U 1254-690, and 4U 1636-536 showed little or insignificant 
spectral variability, thus we could stack their RGS 1 and 2 spec- 
tra, which were taken during all the observations, and produce 
a single final fluxed spectrum for each LMXB (se e Figs. |2] - 
@). We adopted a similar approach in IPinto et ail (1201 Ol) for 
GS 1826-238. There we filtered out the bursting emission from 
the spectra but also reported that it was negligible compared to 
the persistent emission. Therefore we do not remove the periods 
of bursts in any of our data. 

For Aql X-l, SAX J1808.4-3658, and 4U 1735-444 we 
have only found one very good spectrum. Their other observa- 
tions available in the archive have either too short exposure or 
much lower flux in such a way that they do not significantly 
improve the statistics. For each of these sources we simply fit 
the RGS 1 and 2 count spectra simultaneously. The LMXB GX 
339-4 is an exception. Its spectral continuum and slope strongly 
vary, but the three observations with IDs 0148 and 0605 have su- 
perimposable spectra, as well as the remaining two observations 
(IDs 0204). Therefore we split the spectra in these two groups 
and stack them producing two final fluxed spectra with com- 
parable statistics. In Figs. |2]-|4] with show the Ne, Fe, and O 
absorption edges of the GX 339-4 fluxed spectrum obtained by 
stacking the group of observations with ID 0204. 

In Table Q] we report the RGS exposure times and iden- 
tification number for each observation. We also provide the 
Galactic coordinates, distance d and hydrogen column density 
A^h for each source. Most of the distances have been taken 
from [Galloway et al.l (120081) . which estimated them through 
the luminosity peak during the LMXBs burst. The A^h ranges 
refer to the spread between the values measured by the 
Leiden/ Argentine/Bonn Survey of G alactic Hi and the D ickey 
& Lockman Hi in the Ga laxy, see iKalberla et ail d2005l) and 
iDickev & Lockmanl d!990h . The Nu values as provided by these 
two surveys differ by up to 30% and we will determine them 
from the X-ray spectrum. In the RGS fits we rebin the spectra by 
a factor of two, i.e. about 1/3 FWHM (the first order RGS spec- 
tra provide a resolution of 0.06-0.07 A). This gives the optimal 
binning for the RGS and a bin size of about 0.02 A. 

3. Spectral features 

In Figs. 12-14] we show the neutral absorption edges of neon, 
iron, and oxygen for the nine LMXBs, which are the strongest 
spectral features. SAX J1808.4-3658 also shows a prominent 
Ni ls-2p edge at 31.3 A, but most of the sources have weak ni- 
trogen and magnesium edges, which means that the Mg and N 
abundances are not well constrained. The Ne i edge (see 14.3 A 
in Fig. |2j is shallower than those of Fe and O, but it is interesting 
for the presence of additional absorption lines due to mildly ion- 
ized Ne ii-iii and heavily ionized Ne ix gas. For a few sources we 
are even able to detect a weak Nevm line at about 13.7 A. For 
some spectra near the neon edge two weak absorption features at 
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Source 


JO (.a) 


Length (ks) {b) 


l,b M 


d (kpc) {d) 


JV H (10 25 m- 2 ) « 




0060740901 


28.5 








4U 1254-690 


0405510301 
0405510501 


60.6 
60.7 


303.5, -6.4 


15.5 CO 


2.2-2.9 




0105470401 


20.6 








4U 1636-536 


0303250201 
0500350301 
0606070201 


30.9 
31.5 
28.8 


332.9,-4.8 


5.95 CO 


2.6-3.6 


4U 1735^144 


0090340201 


20.6 


346.1,-7.0 


6.5 CO 


2.6-3.0 


Aql X-l 


0406700201 


52.6 


35.7,-4.1 


3.9 CO 


2.8-3.4 


GS 1826-238 


0150390101 
0150390301 


106.3 
91.5 


9.3,-6.1 


6.7 CO 


1.7-1.9 




0148220201 


20.3 










0148220301 


13.7 








GX 339^1 


0204730201 
0204730301 
0605610201 


13.3 
13.4 
32.9 


338.9, -4.3 


8 <s) 


3.7-5.3 


GX 9+9 


0090340101 
0090340601 


10.6 
21.7 


8.5,9.0 


7.5 C" 


2.0-2.1 


SAX J1808.4-3658 


0560180601 


63.7 


335.4,-8.1 


2.77 CO 


1.1-1.3 




0084020401 


21.5 








Ser X-l 


0084020501 
0084020601 


21.7 
21.8 


36.1,4.8 


7.7 CO 


4.0-4.7 



(a,b) identification nu mber and duration o f the observation. < - c - d -^ Source Galactic coordinates in degrees, distan ce from the Sun , and hydrogen 
column density (see iKalberla et al1l2005l and iDickey & Lockman 1990). CO Distance estimates are taken from iGallowav et al] |2008). <s) For 
GX 339-4 we adopt the distance as s ugge sted b y Zdziars ki et al] d2004l) . m The distance for GX9+9 is the mean value between those ones 
provided bv lChristian & Swankl ill 9971) and lZdziarski et al] (12004^ 
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Fig. 1. Map of the X-ray sources as projected on the Galactic 
plane. The Sun is assumed to be 8 kpc far away from the Galactic 
Center (GC). 

15.0 A and 14.2 A are visible that are presumably due to Fe xvn 
and Fe xvm in the hot gas phase. 

The iron L2 and L3 edges are located at 17.15 and 17.5 A 
(see Fig. [3}. The sources show clear differences in the depth of 



these edges. We also provide the position of the Ovii/3 (ls-3p) 
and O vine (ls-2p) absorption lines at 18.6 and 19.0 A, respec- 
tively. These two lines also trace the hot ionized gas. Fig. |4] 
shows the most interesting part of the spectrum: the oxygen 
K edge. The Oi ls-2p line at about 23.5 A is the strongest 
ISM spectral feature, the neutral oxygen is also responsible for 
the jump in spectrum between 22.5 and 23.2 A. Additional but 
not less important ls-2p absorption lines are produced by On 
(23.35 A), Om (23. 1 A), Ovn (21.6 A), and occasionally Ovi 
(22.0 A). Dust oxygen compounds affect mainly the spectral cur- 
vature between 22.7 - 23.0 A (see e.g. lPinto et al]|2010l and ref- 
erences therein). An easy way to check for strong amounts of 
dust is to compare the depth of the edge with that of the line. 
Dust usually is responsible for a deep edge, while the gas also 
provides strong lines. For instance, in Ser X-1 the O K edge is 
deeper than the ls-2p line, while in SAX J1808.4-3658 the ab- 
sorption line is clearly deeper than the edge. This means that 
along the LOS towards Ser X-1 the amount of oxygen found in 
dust compounds is expected to be larger. We will confirm this in 
Sect. [5] Dust is also expected to provide most of iron (see e.g. 
|Jenkinsl l2009). which means that the bulk of the absorption be- 
tween 17-17.5 A is produced by iron dust compounds (see also 
Costantin i et al.l 120121) . Dust is also responsible for most of the 
magnesium edge at about 9.5 A. 
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Fig. 2. Data and best-fitting models: the Ne K edge. 



4. Spectral modeling of the ISM 

The advantage of using LMXBs for the analysis of the ISM lies 
in their high-statistics and simple spectrum. All the sources in 
our sample provided us with well exposed spectra, which are 
easy to model by adopting a simple continuum. In most of the 
cases one power-law (PL) continuum gave an acceptable and 
quickly-converging fit. When a single PL did not provide a good 
fit then we were able to obtain an acceptable fit by simply adding 
another PL. The continuum of LMXBs usually originates from 
their accretion disk and corona in the form of blackbody and 
comptonized emission, respectively. However, we are using only 
the soft part of the X-ray spectrum which is relevant for the 
ISM absorption and is not broad enough to provide accurate con- 
straints on the X-ray emission components. Moreover, this work 
does not focus on the nature of the continuum intrinsic to the 
sources. Thus, we prefer to use PL components rather than phys- 
ical ones like blackbody (BB) or comptonization (COMT). In 
principle the choice of the continuum component may affect the 
estimate of the hydrogen column density, while the other ionic 
columns are well constrained as they are estimated through nar- 



row absorption lines which do not depend on the continuum. 
Thus, we also tested both PL and a combination of BB / COMT 
in modeling our RGS spectra (bb and comt models in SPEX). 
The best fit values of Nh were consistent within the errors, which 
validates our choice of a simple PL continuum. 

We have built an em pirical model for the ISM similar to that 
one successfully used bv lPinto et alj d2012a[) but with a few use- 
ful changes. Essentially, they used three components modeled 
with the slab model of SPEX to fit the three main phases of the 
interstellar gas. The slab model gives the transmission through a 
layer of gas with arbitrary ionic column densities. Here we also 
use two slab components to model the warm and hot phases of 
the ISM, but we prefer to model the cold neutral gas with the 
hot component in SPEX. The hot model calculates the transmis- 
sion of a collisionally-ionized equilibrium plasma. For a given 
temperature and set of abundances, the model calculates the ion- 
ization balance and then determines all the ionic column densi- 
ties by scaling to the prescribed total hydrogen column density. 
At low temperatures this model mimic s very well the neut ral in- 
terstellar gas (see SPEX manual and Kaastra et al. 2009). Free 
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Fig. 3. Data and best-fitting models: the Fe L2 and L3 edge. 



parameters in the hot model are the hydrogen column density 
Nu, the temperature T, the velocity dispersion cr v , the systematic 
velocity v, and the abundances. However, in X-rays the spectral 
resolution is not high enough to constrain the velocities with suf- 
ficient accuracy, thus we assume that the interstellar gas is at rest. 
We also assume a temperature of 0.5 e V (about 5800 K), the low- 
est value available in the SPEX hot model, because in this way 
the gas is mostly neutral. In X-rays with the current satellites 
it is not possible to resolve the narrow lines of the interstellar 
medium. A spectral fit can easily provide a wrong value of veloc- 
ity dispersion cr v and consequently of column density because 
they are degenerate. Therefore it is more appropriate adopting 
a physically reasonable value for the cr v and just fitting the col- 
umn density of each ion. We adopt a nominal value of 10 km s 
for the velocity dispersion of the cold (O i) and wa rm (Oii-v) 
gas co mponents, which is similar to that found by iPinto et aD 
(1201 2al) using a higher-resolution UV spectrum. For the hot gas 
(Ovi and higher ionization states) we adopt cr v =100kms _1 , 
which is an average value between th ose suggested in th e lit- 
erature (see for instance the work by lYao & Wangl 120051 and 



lYao et al.l 120091 on high-resolution Chandra X-ray spectra of 
several LMXBs). We took into account absorption by interstel- 
lar dust with the SPEX amol component. The amol model cal- 
culates the transmission of various mole cules, for details see 
IPinto et all d2010l) . ICostantini etall (l2012h and the SPEX man- 
ual (Sect. 3.3). The model currently takes into account the mod- 
ified edge and line structure around the O and Si K-edge, and 
the Fe K / L-edges, using measured cross sections of various 
compounds taken from the literature. Relevant and abundant sili- 
cates, ices and organic molecules are present in the our database. 
In fitting the molecular models we limit the range of column 
densities to physical values. For instance, we took care that the 
column density of each weakly-abundant element did not ex- 
ceed the proto-Solar value predicted by the best-fitting A^h- This 
is important especially for Ca and Al (contained respectively in 
andradite and hercinite), whose abundances are much lower than 
oxygen. 
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Fig. 4. Data and best-fitting models: the O K edge. 



4.1. The fitting procedure 

The first step is determining a preliminary model for the whole 
spectrum. This requires to fit the spectrum of each source with 
a simple model consisting of one (or more) power-law compo- 
nents) absorbed by neutral gas (hot model in SPEX). As already 
mentioned, one or two power-law components are sufficient to 
reproduce all the spectral continua. This fit provides prelimi- 
nary values for the normalization and slope of the power-law 
as well as the hydrogen column density and the abundances of 
the neutral gas. The second step is determining the contributions 
of each absorbing component. Some compounds and ions pro- 
duce relevant absorption features only near certain absorption 
edges. For instance, ions such as On-vi and molecules like wa- 
ter ice, andradite and pyroxene mostly affect the spectral range 
near the O K-edge. Therefore, we first fit each edge with the 
preliminary model and add in sequence ions and, in a second in- 
stance, molecular compounds looking for significant changes in 
the chi 2 . In the case of the oxygen K-edge, significant changes 
in the x 2 are given by ions like O ii-iv, O vii-vm, and compounds 
like water ice, pyroxene, magnetite, and hematite. Once we have 



measured the Ax 2 of the absorbing components in each edge, we 
jointly fit the edges or the spectral ranges that share transitions 
from same ions or compounds. The wavelength interval between 
15 and 24 A contains the strong Fe and O edges, and the bulk 
of the dust signatures. The fit of this spectral region provides a 
starting model for the final fit extended to the whole spectrum. 

4.2. Step-by-step analysis for GX 339-4 

In this section we show the detailed analysis of the two GX 339- 
4 stacked spectra between 15 and 24 A. This will also indicate 
the relevance of each absorber in the final model. We choose 
to show this source because of its highest statistics and large 
column density that maximize the significance of the spectral 
features. Moreover, the results obtained for GX 339-4 are sim- 
ilar to those of most sources in our sample. We first jointly 
fit the Fe L and O K edges (15 - 24 A) with a simple model 
consisting of a power-law absorbed by neutral gas. Free pa- 
rameters are the normalization and slope of the power-law and 
the Fe and O abundances of the neutral gas. This fit provides 
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Table 2. GX 339-4 spectral fits within 15 
A X 2 . 



24 A and relative 



Parameter 


A* 2 


Parameter 


A* 2 


Parameter 


A* 2 


On 


14 


Metallic Fe 


101 


CO 


1 


Om 


163 


H 2 O ice 


168 


FeO(OH) 


1 


Oiv 


20 


MgSiO a 


4<«) 


N 2 





Ov 


1 


Fe 3 4 


60 


Fe 2 3 





Ovi 


1 


Ca 3 Fe 2 Si 3 12 


3 


FeS0 4 





Ovn 


366 


FeAl 2 4 


2 






vm 


35 


Mgi. 6 Fe .4 SiQ 4 


2 







The Ax 2 values refer to the sequential changes obtained by adding the 
parameters from top to bottom and from left to right: O n, O m, O vm, 
Metallic Fe, FeS0 4 . (n> The A% 2 given by pyroxene MgSi0 3 is 
small because its features are similar to those produced by H 2 O ice 
(see also Fig.|7J. 



xl - 1692/454 = 3.73, which is clearly poor and unacceptable 
(see also the blue line in Fig. [5). A different continuum does not 
produce any change. Therefore, we add ions in sequence from 
O ii to O vm, see Table|2] This (neutral + ionized) gas model with 
all these ions gives xl = 1092/447 = 2.44 and a big improve- 
ment to the fit. However, the fit is not yet acceptable and there 
are residuals at 17.5 and 22.9 A (see red line). Moreover, the 
pure gas model does not reproduce the ratio between the absorp- 
tion lines and edges for both iron and oxygen. Dust is definitely 
needed. We complete the ISM model by adding in sequence dif- 
ferent molecules with the order as suggested by preliminary fits 
to each edge, which is mentioned above: metallic iron, water ice, 
pyroxene silicates, up to iron sulfate (see also Table|2}. The final 
xl drops down to about 1 .7 as shown by the final model (purple 
line in Fig. [5j. Although this fit is not yet formally acceptable, 
the final xl (significantly higher than one) is mainly due to the 
high statistics. Similar;^ 2 values are obtained for featureless re- 
gions at similar statistical quality and indicate the final system- 
atic uncertainties of the instrument model. In Fig. [6] we show the 
transmission of the ISM best-fitting model near the O i edge of 
GX 339-4. We also report the relative contributions provided by 
O i (red solid line), O n (purple dotted line), O iii-iv (black dashed 
line), as well as by specific compounds. Molecular compounds 
dominate the spectral range within 22.8 - 23.2 A and affect the 
ratio between the depths of the O i edge and resonance line. 

Once we obtain a good fit for the 1 5 - 24 A spectral range, we 
fit the model to the entire spectrum. Additional free parameters 
in the model are the neon and magnesium abundances of the 
cold gas (hot model in SPEX) and the column densities of Ne n- 
iii, Ne vm-x, Mg xi- xii, and Fe xvn-xvm (which are provided by 
the slab model). 

4.3. Results 

This procedure used for GX 339-4 is identically applied to the 
other eight targets and the results of our best-fitting models ob- 
tained on all the X-ray sources are shown in Table [3] and plotted 
in Figs. [2}j4] The model fits very well both the absorption edges 
and lines. The reduced^ 2 for the nine sources are between 1.4 
and 1.7. We show the absolute abundances and hydrogen col- 
umn density for the cold gas, the column densities for all ions 
and molecules that we were able to detect. In the cases of no 
detection we report the 2 cr upper limits. We also report the for- 
mula for each compound and the total amount of iron and oxy- 
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Fig. 5. Spectral fits of the Fe L and O K edges for GX 339-4 (for 
clarity only one stacked spectrum is shown). For more detail see 
Sect. l4~2l and Tablel2l 
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Fig. 6. Transmission of the ISM near the O K-edge in the spec- 
trum of GX 339-4. The model refers to the best-fitting model as 
given in Table[3] The solid red and black lines provide the trans- 
mission of the neutral gas and the entire ISM, respectively (see 
also Sect.l4~2l. 



gen found in dust. We will discuss the results in Sect. Briefly, 
we have found a general agreement in the spectra of the dif- 
ferent sources. Ions like O v-vi and molecular compounds like 
hematite (Fe2C>3), and CO are hard to detect. Instead, O vn-vm, 
H2O ices, and compounds of Ca (andradite), and Al (Hercinite) 
are detected along any LOS. The cold neutral gas provides on 
average Oi ~ 1 - 3 x 10 22 rrT 2 , which means that most oxygen 
is found in the neutral phase. We have summed all the contri- 
butions to the molecular oxygen as well as all the ionic column 
densities of the mildly ionized (0 11- v) and heavily ionized (O vi- 
vm) gas, and compare them in Table [4] As we have mentioned 
above, most of the oxygen is provided by the neutral gas phase. 
Dust also seems to be ubiquitous. The warm and hot ionized gas 
phases generally contribute less to oxygen, but there are some 
cases in which they constitute significant portions of the total 
column density. SAX J1808. 4-3658 is an exceptional case with 
detectable O vi and prominent O vn absorption lines (see also 
Fig. |4l bottom-right panel). 



Table 3. RGS fits with the complete ISM model. The dust and the three gas phases are represented by the four main blocks. 



oo 



Parameter 


4U 1254-690 


4U 1636-536 


4U 1735-444 


Aql X-l 


GS 1826-238 


GX 339^1 


GX 9+9 


SAX J1808.4-3658 


SerX-1 




N H (10 25 m" 2 ) 


2.69 ± 0.03 


3.58 ± 0.07 


3.28 ±0.08 


5.21 ± 0.05 


4.14 ±0.06 


5.1 ±0.4 


2.15 ±0.05 


1.40 ±0.03 


5.0 ± 0.3 




0/H (,,) 


0.92 ± 0.04 


0.85 ±0.10 


0.93 ± 0.03 


0.88 ± 0.02 


0.92 ± 0.04 


1.06 ±0.05 


0.99 ± 0.03 


0.91 ±0.03 


0.9 ±0.1 




Ne/H (fl) 


0.91 ± 0.03 


1.31 ±0.05 


1.00 ± 0.05 


0.95 ± 0.07 


1.23 ±0.05 


1.05 ±0.05 


1.35 ±0.07 


1.37 ±0.05 


1.20 ±0.05 




Mg/H<°> 


1.7 ±03 


<0.8 


2-0 ±?n 

l.U 


U.a 


1-4 ±™ 

0.6 


1.2 ± j* 

O.o 


0.4 + 0-5 

0.3 


<0.2 


U. 1 




Fe/H ( "» 


0.16 ± °f 

(J. 1 


0.18 ±0.04 


0.25 ±0.10 


0.11 ±0.05 


0.22 ± 0.06 


0.4 ±0.1 


0.19 ±0.05 


<0.1 


0.50 ± 0.05 




On« 


0.33 ±0.15 


1.4 ±0.7 


0.6 ± 0.2 


0.2 ±0.1 


0.6 ±0.1 


1.4 ±0.5 


0.3 ±0.1 


1.67 ±0.05 


2.0 ± 0.2 




Om 9) 


0.08 ±°™ 

0.04 


<0.1 


0.2 ±0.1 


<0.3 


<0.3 


<0.8 


<0.05 


0.34 ± 0.05 


<0.4 




Oiv<« 


0.5 ± 0.2 


0.3 ± 0.2 


0.4 ± 0.2 


0.4 ±0.1 


0.6 ± 0.2 


0.5 ±0.1 


0.4 ± 0.2 


0.31 ±0.03 


0.7 ± 0.5 


CP 


Ov< 6) 


<0.04 


<0.02 


<0.1 


<0.03 


<0.06 


<0.01 


< 0.002 


0.006 ± 0.002 


0.3 ± 0.2 


Neil*' 


0.11 ±0.06 


1.2 ±0.3 


0.4 ± 0.2 


0.7 ± 0.2 


1.4 ±0.3 


1.2 ±0.5 


0.4 ± 0.2 


1.3 ±0.1 


0.3 ± 0.2 


5 



Nem' 6 ' 


0.19 ±0.08 


0.9 ± 0.2 


0.5 ± 0.2 


0.6 ± 0.2 


0.6 ± 0.3 


0.9 ± 0.2 


0.2 ±0.1 


1.1 ±0.1 


0.8 ± 0.4 


re 


Ovi IW 


<0.03 


<0.06 


<0.01 


<0.06 


<0.07 


<0.04 


< 0.004 


0.02 ± 0.01 


<0.05 


& 


Ovn (i) 


0.10 ± 0.05 


0.20 ± 0.07 


0.2 ± 0.1 


0.15 ±0.06 


0.15 ± 0.06 


0.4 ± 0.1 


0.04 ± 0.02 


2.2 ± 0.5 


0.24 ± 0.05 


ISM c( 


Ovm"" 


0.09 ± 0.04 


0.19 ± 0.05 


0.3 ± 0.1 


0.10 ±0.04 


0.10 ± 0.03 


0.3 ± 0.1 


< 0.005 


0.11 ±0.02 


0.20 ± 0.06 


Nevm**' 


<0.03 


0.07 ± 0.04 


<0.05 


<0.01 


<0.01 


0.06 ± 0.04 


< 0.004 


< 0.003 


<0.03 


1 


Nerx® 


0.05 ± 0.03 


0.09 ± 0.04 


0.2 ±0.1 


0.07 ± 0.05 


0.06 ± 0.02 


0.16 ±0.04 


0.03 ± 0.02 


0.05 ± 0.02 


0.14 ±0.05 





Nex (i) 


<0.04 


<0.03 


<0.02 


<0.1 


<0.02 


<0.02 


< 0.006 


< 0.008 


<0.04 


§ 


Mgxi (i ° 


0.16 ±0.05 


0.08 ± 0.04 


0.4 ± 0.3 


<0.05 


<0.1 


0.10 + 0.05 


0.06 ± 0.03 


0.18 ±0.09 


< 0.15 




Mgxn (i) 


<0.01 


<0.01 


<0.01 


<0.02 


<0.02 


<0.01 


< 0.005 


<0.01 


<0.03 


r 


Fe xvii (6) 


<0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.008 ± 0.004 


0.02 ± 0.01 


<0.01 


|h X-ra 


Fe xvm 


<0.02 


<0.01 


<0.01 


0.02 ± 0.01 


<0.01 


<0.01 


< 0.003 


<0.04 


<0.01 


Water ice H 2 O (b) 


2 ± 1 


0.9 ± 0.5 


3 ± 2 


3 ± 1 


3 ± 2 


4.4 ± 2.5 


2 ± 1 


0.3 ±" 

— Q J 


< 6.6 


y spectres 


Andradite Ca 3 Fe, Si 3 12 (b) 


0.02 ± 0.01 


0.02 ± 0.01 


0.02 ± 0.01 


0.03 ± 0.01 


0.03 ± 0.01 


0.03 ± 0.01 


0.012 + 0.004 


0.008± 0.004 


0.03 ± 0.01 


HerciniteFeAl 2 4 {b) 


0.04 ± 0.02 


0.05 ± 0.02 


0.04 ± 0.02 


< 0.025 


0.06 ± 0.02 


0.07 ± 0.02 


0.03 ± 0.01 


0.015± 0.007 


0.07 ± 0.02 


Lepidocrocite FeO(OH) {b) 


0.11 ± 0.06 


0.4 ± 0.2 


< 0.2 


0.6 ± 0.3 


0.4 ± 0.2 


0.03 ± 0.02 


0.2 ± 0.1 


0.15 ± 0.08 


0.15 ± 0.10 


cop 


Pyroxene Mg Si 3 m 


< 1.2 


1.2 ± 0.5 


< 1.1 


< 1.2 


< 1.1 


0.3 ±°-' 

0.1 


< 1.0 


< 0.4 


1.7 ± 0.5 


•~< 



Carbon monoxide CO 


< 0.11 


0.3 ± 0.2 


< 0.2 


< 0.1 


< 0.3 


0.07 ± 0.05 


< 0.05 


0.6 ± 0.1 


< 0.2 


— 1 
r 


Metallic Fe <b) 




°- 3± ai 


0.3 ±g; 2 


<0.7 


< 1.2 


0.6 ±0.1 


°- 3± Si 


0.13 ±0.05 


<0.7 




Magnetite Fe 3 4 (6) 


0.03 ± 0.02 


0.2 ±0.1 


<0.04 


0.3 ± 0.2 


0.3 ±0.1 


0.2 ±0.1 


<0.02 


< 0.005 


0.2 ±0.1 




Olivine Mg L6 Fe 4 Si0 4 m 


<0.04 


0.06 ± 0.03 


<0.1 


<0.04 


<0.08 


0.07 ± 0.03 


0.2 ±0.1 


<0.03 


<0.1 




Laughing gas N 2 


0.4 ± 0.3 


<0.04 


<0.1 


<0.03 


<0.2 


<0.03 


<0.03 


0.03 ± 0.02 


<0.04 




Hematite Fe 2 3 {b) 


<0.04 


<0.01 


<0.1 


<0.02 


<0.2 


< 0.003 


<0.02 


<0.01 


<0.02 




Iron sulfate Fe S0 4 ib) 


<0.04 


<0.01 


<0.05 


<0.03 


<0.01 


< 0.004 


<0.04 


<0.03 


<0.02 




Total oxygen in dust {b> 


3 + 1 


7±2 


4± 1 


6± 1 


6± 1 


7 ±2 


4± 1 


1.5 ±0.5 


7 ±2 




Total iron in dust ™ 


0.6 ± 0.2 


1.3 ±0.5 


0.5 ± 0.3 


1.7 ±0.5 


1.5 ±0.6 


1.4 ± 0.5 


0.7 ± 0.2 


0.3 ±0.1 


0.8 ± 0.4 





Abundances ratios are in the linear proto-Solar abundance units of Lodders & Palme (2009). 
Ionic and molecular column densities are in units of 10 2 ' nr 2 (see also Sect.[4j- 
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4.4. Limitations 

There are some factors that may limit and affect the uniqueness 
of our best-fitting models. First of all, the features produced by 
molecules in a certain absorption edge are similar and degen- 
erate. Near the Fe L-edge the RGS absorption lines blend and 
may not give unique solutions. Chandra spectra are provided 
with higher resolution, but they give accurate results only for a 
few sources even brighter than ours. Therefore, in our standard 
modeling we have simultaneously fitted all the molecules and 
all the absorption edges and lines in order to obtain the best- 
fitting mixture and to break as much as possible the degeneracy. 
Unfortunately, for some compounds, such as magnetite and py- 
roxene, we only have the transmission near the oxygen edge, 
while for the other edges we simply use the pure atomic cross- 
section without absorption lines (for details consult the amol 
model in the SPEX manual). Fortunately, these systematic ef- 
fects do not strongly alter the depletion factors, especially for 
oxygen. The Mg K-edge is very weak and it does not make 
any difference whether we model it with atomic or molecular 
cross-sections. Another problem is the degeneracy between the 
absorption features produced by water i ce and pyroxene (see 
Fig. [6] and [7] and ICostantini et al.l 1201 2l Appendix A), which 
are the molecules that best fit the O K-edge. This is also shown 
by the large uncertainties on their column densities in Table [3] 
Therefore, in order to test the robustness of our results and to es- 
timate the systematic effects, we have performed additional fits 
with different combinations of molecules by excluding some of 
them. One model includes only the molecules that do not lack 
any relevant cross-section, such as water ice, Metallic Fe, CO, 
and Hematite (all the other molecules such as pyroxene and mag- 
netite are excluded). In another test we replace water ice with py- 
roxene in order to measure systematic changes in the Mg abun- 
dance of the cold gas. Several tests like these provided us with 
systematic changes in the column densities of the molecules and 
the ions as well as in the abundances of the cold gas. We have 
summed all the systematic and statistic errors and reported them 
together with the best-fitting values in Table [3] Interestingly, the 
total amount of oxygen and iron in dust does not strongly depend 
on the chosen dust mixture. Therefore, their depletion factors are 
not highly effected, but the detailed chemical structure is. 

On another hand, the molecular database is not complete. We 
have tested all the molecules available in the SPEX database, 
which are about 40 (carbon oxides, ices, iron and magnesium 
silicates, and more complex molecules). Although these com- 
pounds are among the most abundant in the ISM, we could still 
miss important co ntribution f rom other species like forsterite 
(Mg 2 Si 4 , see e.g. Uonesll2000l) . Our ISM team at SRON is cur- 
rently taking synchrotron and electroscopic measurements of the 
cross-sections near the O, Mg, and Si K-edges, and both the Fe 
K/L edges of several silicate compounds with a high spectral 
resolution too. These laboratory data will provide accurate esti- 
mates of molecular column densities. 

The results obtained on the various absorption edges have 
different weights. The oxygen edge is the deepest and the col- 
umn densities of oxygen compounds are constrained better than 
for iron. The Mg edge is very weak for hydrogen column den- 
sities like those of our sources and the estimates of magne- 
sium abundances and depletion factors are much more uncer- 
tain. This does not apply to neon as it is supposed to be only 
in a gaseous phase. This work mostly focuses on measurements 
of interstellar oxygen, iron, and neon column densities. A deep 
analysis of magnesium and silicon would require the use of 
both the XMM-Newton RGS and the Chandra HETG gratings 




N 



(10 21 rrT 2 ) 



Fig. 7. Water ice / pyroxene anti-correlation. The pyroxene col- 
umn densities refer to the spectral fits of the Fe L and O K 
edges for GX 339-4 obtained with fixed amount of water ice. 
The red values show the x 1 values of each fit. For more detail 
see Sect.l4~2landl4~4l 



(which have higher effective area and resolution than RGS at 
short wavelength) in combination with sources with higher hy- 
drogen column densities. 

Our estimates of column densities and abundances provide 
interstellar properties as integrated along the LOS and on a large 
scale. This is reasonable for the diffuse gas, which is almost uni- 
formly distributed between 3-15 kpc from the Galactic Center 
and whose warm phase may reach a few kpc height from the 
Galactic place (see e.g. Ferriere 2001). This scale-height is in- 
deed comparable to the average altitude of our sources (see 
Table [1). However, this assumption may not be correct for dust 
and molecules, which follow the Galactic spirals and are usually 
found at lower altitudes. This means that dust sampling towards 
our targets may be closer to the Sun rather than in the middle 
of the path. In the cold phase, iron and magnesium are strongly 
depleted into dust. Therefore, the average location of the Fe and 
Mg absorbers can be found at shorter distances than those refer- 
ring to oxygen and neon. However, in Sect. l5.2l we will show that 
along our low-latitude LOS the different scale-heights of gas and 
dust will not affect the results. 

Another problem in our grating spectra is the presence of 
several bad pixels or dead columns that affect the estimates of 
some parameters. One of the most relevant is located at 22.75 A, 
just near the Oiv resonance line (see Figs. |4] and |6}. The Oiv 
column densities of our targets are all determined just through 
this line and their values may be inaccurate. New observations 
taken with the multi-pointing RGS mode will surely provide a 
solution to this problem. 



5. Discussion 

Our analysis shows that the ISM has a multi-phase structure 
characterized by cold neutral gas (and dust), mildly ionized gas 
and heavily ionized gas. This complex structure is found to- 
wards any of the studied source s (see Table [3) and is consistent 
with our previous work (see e .g. Kaast raet al.l l2009: Pinto et all 
2010. 1201 2at ICostantini et al.ll2012l) . Signatures of dust are al- 
ways present with main contributions by ices and silicates. The 
oxygen edge and lines provide the most accurate results due to 
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their large depth, but the results are even more robust because 
of the simultaneous modeling of all the absorption edges. This 
was crucial especially for fitting dust and molecules which af- 
fect both O, Fe, and Mg edges. An interesting and simple check 
of different amounts of oxygen in dust is provided by the com- 
parison of the OK edge in Fig. [4] As we anticipated in Sect. [3] 
the O i ls-2p absorption line is deeper than the oxygen edge for 
SAX J1808.4-3658, while for Ser X-l and GX 339^ they are 
comparable. We attribute this to oxygen depleted from the gas 
phase into dust grains, which is confirmed by the measurement 
of a larger amount of oxygen in solids in the LOS towards the 
two latter sources (see Table |3). 

5.1. The nature of the gas phases 

In our standard model we have used a collisionally-ionized 
model in SPEX because at low temperatures this mimics 
well the cold interstellar gas. However, at these tempera- 
tures photo-ionization could also provide a significant contri- 
bution. Therefore, we tried an alternative model in which a 
photo-ionized component {xabs in SPEX) substitutes the cold 
collisionally-ionized hot component. The fits are equivalent, the 
gas is almost neutral with a very low ionization parameter 
Abundances and^- 2 are consistent with our standard model. The 
main difference is that the fit with the photo-ionized model takes 
more time to converge and the absolute abundance s, i.e. rela- 
tive to hydrogen, have larger uncertainties (see also iPinto et all 
I2012al) . Therefore, we prefer to keep our standard model even if 
photo-ionization is not ruled out. 

We have also considered an alternative, physical model for 
the warm mildly-ionized gas that produces most of the On-iv 
and Neii-m absorption lines. Instead of a slab absorber we have 
tested a photo-ioniz e d xabs component, which was successfully 
used by IPinto et aTl d2012al) in the deep UV and X-ray spectra 
of the AGN Mrk 509. The fit was poor because a single photo- 
ionized component is not able to fit all the large column den- 
sities as measured by our empirical standard model and shown 
in Table[3] We need a low-ionization component which provides 
the bulk of O n and Ne n together with an intermediate-ionization 
component for the other ions. Interestingly, a good fit is ob- 
tained when using the same abundance pattern as provided by 
the cold neutral gas, while proto-Solar abundances give a bad fit. 
However, a thorough study of the warm gas requires longer ex- 
posures and complementary UV data, which is beyond the scope 
of this paper. 

The heavily ionized gas is usually found to be in collisional 
equilibrium. However, it is thought that this hot gas in the Galaxy 
is not characterized by a single phase, but it is constituted of 
two main phases. Most Ovn-vm should arise from a very hot 
(~ 2 x 10 6 K) collisionally-ionized phase, while the bulk of 
Ovi is thought to be embedded in a cooling phase with tem- 
peratures 1 - 5 x 10 5 K (see e.g. Richter 2006). We have tested 
both single-phase and two-phase models by substituting the slab 
component responsible for the O vi-viii absorption with one and 
two collisionally-ionized hot components in SPEX. The fits are 
comparable due to the weakness of their absorption lines. We 
would need UV O vi data in order to put strong constraints, but 
most of our targets have not been observed with COS/HST or 
FUSE. However, Table [3] shows a spread larger than 10 in the 
Ovu / Ovm column ratio, which for a single gas component 
would provide a very large scatter in temperature between the 
nine LOS (from 1.2 to 2.6 million K). We think that combina- 
tions of two stable, but different hot (cooling and coronal) phases 
is a more likely description. 



5.2. Is the ISM chemically homogeneous ? 

On small scales the ISM is not expected to be highly homoge- 
neous because there are significant physical and chemical dif- 
ferences between the interstellar environments like H n regions, 
PNe nebulae, dark clouds, etc. However, the diffuse medium in 
our Galaxy may be homogeneous on large scales. Some chem- 
ical and physical properties may be possibly consistent when 
comparing distant regions which underwent a similar evolution. 
This is possible for dust as most o f it is expected to grow in th e 
ISM rather than in stars (see e.g. Mattsson & Andersenll2012l) . 
It is thus worth to compare its chemical properties as integrated 
along the different LOS. First of all, we notice that the molecular 
composition is similar in any direction. Aluminates and calcium 
silicates are found to provide most of Al and Ca, which are in- 
deed expected to be highly depleted into dust. Water ice may 
be present, but in an uncertain amount du e to the degen eracy 
between its features and those of pyroxene. Uenkinsl d2009) sug- 
gested water ice as a possible principal resource of oxygen in 
dust, but he also showed that at low-intermediate Nh it is hard to 
detect. Cross-sections of these and of additional compounds are 
currently being measured with a resolution higher than the cur- 
rent one and may provide more accurate results (see Sect. I4.4I >. 
To understand the role of dust, we have summed all the contri- 
butions to oxygen from solids and molecules and compared it 
with the hydrogen column density that we have obtained in our 
best fits (see Fig. [8}. We show the trends of both the total column 
density of dusty oxygen and its fraction with respect to the total 
amount of neutral oxygen (gas + dust) as function of the Nh. The 
oxygen dust column correlates well with the hydrogen column 
density (slope = (1.3 + 0.2) x 10~ 4 ) and the fraction of oxygen 
bound in dust is constant (19 + 5%) along any LOS. These re- 
sults suggest that on average the cold ISM phase is chemically 
homogeneous. We notice that the total amount of oxygen in dust 
is solidly measured as it does not strongly depend on the dust 
mixture. For instance, we tested different molecular combina- 
tions by removing some of the main compounds from the fit, 
but in the end the sum of all the oxygen molecular column den- 
sities provide the same results in the fits. The same applies for 
iron, whose dust fraction does not depend on the Nh and the 
LOS, but it has a larger uncertainty and spreads between 65- 
90%. This occurs because the iron L edge is shallower than the 
oxygen K edge. As we mentioned in Sect. 14.41 the scale-height 
of the interstellar dust is smaller than for the gas, thus, due to the 
altitudes of the targets, the dust is expected at shorter distances 
from the Sun. This means that the analysis of the dust in these 
different LOS probes properties of the ISM as measured on a 
large scale, but probably smaller than with the gas. As additional 
check we have computed the dust-to-gas ratio for our sources as 
the ratio between the total amount of oxygen in dust and the hy- 
drogen column density. The dust-to-gas ratio is constant (slope 
= 0.0 + 0.2 x 10~ 4 ) at any value of Galactic latitude as shown in 
Fig. [9] This result and the trend seen in Fig. [8] suggest that the 
distribution of the dust should not deviate from the gas. This may 
be due to the low latitudes of our sources, which was also one 
of the reasons for us to choose these targets as representatives of 
the Galactic (inner) disk. 

5.3. ISM column densities 

Our empirical model provides us with estimates of column den- 
sities for several molecular and ionic species (see Table [3). In 
some cases we could only obtain upper limits, but we prefer to 
show them as they may be useful for future work. In order to fur- 
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Fig. 8. Amount of O i in dust and its fraction with respect to the 
total neutral oxygen (gas + dust) versus the hydrogen column 
density. 
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Fig. 9. Ratio between the Amount of O i in dust and the hydrogen 
column density versus the absolute value of the Galactic latitude. 

ther check the robustness of our model we compare our column 
density estimates with those found in the literature. In most cases 
our best-fit values of Nh agree with those measured at longer 
wavelength (see Table Q}. On average our values are slightly 
higher that the 21cm Hi estimates probably because the cur- 
vature as measured in the X-ray above 25 A takes into account 
the contributions fro m H i and H2. Our values are consistent with 
other papers (see e. g. Xiang et al. 2005; Caballero-Garcfa et al. 
2009; Bh attacharvv a & Stro hmaverll2007l iPatruno et all 12009k 
Our column density estimate for G S 1826-238 agrees with what 
we found with a different model in Pinto et al . (2010), i.e. N H = 
4.14 ± 0.07 x 10 25 irT 2 . All the neutral column densities that 
we have measured for 4U 1636-536, 4U 17 35-444, GX 339- 
4, GX 9+9, and Ser X-l are consistent with lJuett et al.1 (120041 
2006). We found a disagreement between our average values 
of Oii/Oi and Om/Oi column densit y ratios, ~0.06 and ~0.02 
respectively, with those obtained by lJuett et aD (120041) of about 
0. 1 . This may be due to the fact that their model does not take 
into account any contribution by dust or molecules near 23 A 
and thus may overestimate the amount of ionized oxygen. They 
obtained acceptable fits without including dust in their models, 
but this was due to their short exposures and to the sensitivity 



of Chandra which is smaller than XMM-Newton near the O K- 
edge. Our ionic column densities agree very well with those es- 
timated bv lMilleret al.l(l2004l) for GX 339-4. According to our 
best fit models, O iv appears to have column densities larger than 
O in. The most obvious cause may be the presence of a bad pixel 
exactly near the O iv resonance line as mentioned in Sect. 14.41 
However, it may also be that we have slightly over-estimated the 
absorption of molecules and dust near at 21 A at expensive of 
the O m amount. Recently, the XMM-Newton/RGS has been pro- 
vided with an optimal observation scheme, the multi-pointing 
mode, which allows snapshots of the same source with shifts 
in the dispersion direction. These shifts place the bad pixels in 
different location and allows to reconstruct a completely clean 
spectrum. This procedure has been su ccessfully used in large 
campaign on AGN Mrk 509 (see e.g. iKaastra etaD 120 1 lh and 
provided our team with a bright and clean spectrum. Once ap- 
plied to sources with high N m H like LMXB, this procedure 
will reveal the composition of the ISM in detail. We point out 
that many of the high-resolution X-ray (RGS) spectra of these 
sources have been analyzed for the first time by us. 

5.4. ISM abundances 

We have also estimated the total a bundances of O, N e , Mg , and 
Fe for the cold phase in units of lLodders & Palmd (120091) by 
summing, for each element, the contributions from neutral gas, 
dust, and molecules (see Table|5]). In most cases all the elements 
appear to be over-abundant with respect to the proto-Solar val- 
ues. This is expected as all LOS cross inner Galactic regions 
where the abundances are higher. It is very difficult to estimate 
the location of the absorbing medium for several reasons. At 
first, we integrate over the LOS and thus our abundances have 
to be considered as large-scale properties of the ISM. Moreover, 
the ISM mass-density distribution within the Galaxy is not well 
known and the LMXB distances are affected by a large uncer- 
tainty which will dominate the error on the ISM absorber loca- 
tion. Therefore, we adopt half of the distance to each LMXB 
as a first-order, average-distance indicator. This means that the 
Galactocentric radius Rq of each ISM absorber will be defined as 
the distance from the Galactic Center (hereafter GC) to the mean 
point of the segment which connects the Sun to each target. We 
show the abundance trends for O, Ne, Mg, and Fe with respect 
to Rq (in kpc) in Fig. [TO] These trends were fitted with a linear 
function and we obtained the following abundance gradients: 

12 + log(0/H) = (8.97 + 0.04) -(0.021 +0.008) xR G , 

12 + log(Fe/H) = (7.76 ±0.14) -(0.034 + 0.026) x R c , 

12 + log (Ne/H) = (8.34 + 0.04) - (0.034 + 0.007) x R c , 

12 + log(Mg/H) = (7.87 + 0.12) -(0.026 + 0.023) x R c . 

The plots show that only the oxygen abundance is clearly anti- 
correlated with the distance from the GC. Neon also shows 
a gradient but with a large scatter probably due to contribu- 
tion i ntrinsic to the so urces (see e.g. lJuett & Chakrabartvll2003l 
and lMadei et al.ll20loh . However, the linear fits show negative 
gradients for all the abundances which is consistent with the 
fact that stellar evolution proceeds faster in the Galactic in- 
ner regions. Interestingly, our abundance gradients agree with 
those estimated with different methods and in different wave- 
length regimes dMaciel & Ouirezalll999t iRolleston etail 2000: 
lEsteban et al] 12005: Ped icelli et al.ll2009h . On average our val- 
ues are slightly smaller than those found in the literature. An 
explanation may be that we overestimated the distance of the 
absorbers by a (small) systematic factor. In fact, the LMXBs are 
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not exactly in the Galactic Plane, but are above or below by some 
hundreds parsec. This means that the distance of each ISM ab- 
sorber, as projected on the Galactic plane, should be smaller than 
what we have adopted. The altitudes of the LMXBs also affect 
the abundance estimates because of the (nega tive) vertical metal- 
licity gradient (see e.g. [Pedicelli et al.ll2009h . The slopes of the 
radial an d vertical gradients are similar (we have compared their 
effects in lPinto et al.ll2012a) . The sources have low latitudes (see 
Tabl^D, on average we obtained b ~ 6.1°, which implies that the 
radial abundance variations are about 10 times larger than the 
vertical ones. This explains why ignoring the vertical gradient 
may only marginally affect our results. We may have underes- 
timated the Mg and Fe abundance gradients also because in the 
cold phase these elements are locked up into dust. In fact, as re- 
ported in Sect. 14.41 dust could be found at distances shorter than 
the gas along our LOS. This means that the average distances of 
the Fe and Mg absorbers may be smaller and the gradients higher 
than those estimated. However, once taken all this into account, 
our estimates of abundance gradients agree even better with the 
literature. 

We can also predict the local proto-Solar interstellar abun- 
dances and check if they agree with those found in the Solar 
system. According to the equation above, we calculate that at 
the Galactoce ntric radius of the Sun, which is ~8.0kpc (see e.g. 
Malkin 2012), the predicted proto-Solar abundances are 

12 + log(0/H) = (8.80 + 0.08), 
12 + log(Fe/H)Q = (7.48 + 0.25), 
12 + log(Ne/H) = (8.07 + 0.07), 
12 + log(Mg/H)Q = (7.66 + 0.22). 

These values agree very well with the abundanc es in the Solar 
System as measured bv lLodders & Palmel (120091) [0/H]= 8.76, 
[Fe/H]= 7.54, [Ne/H]= 7.95, [Mg/H]= 7.62, and further validate 
our method. Interestingly, the interstellar oxygen abundance, as 
measured by previous work, is always lower than the abundance 
estimates for the Solar System, OB stars, H n re gions, planetary 
nebul ae, Cepheids, and Sun-like s tars (see e.g. Stasiriska et alJ 
120121 and iNieva & Prz vbillal [2012b . Instead, our measurements 
and predictions agree very well with the abundances of these 
stellar indicators, which confi rms that the local Gal axy is chem- 
ically homogeneous (see also Stasihska et ai] |2012l) . 

As clearly seen in Fig. [10] our work is optimal to probe ISM 
abundances within 4-8 kpc from the GC, but this range could be 
easily extended by using more sources. This may be an interest- 
ing proxy for a future work extended to more LOS and back- 
ground source types like LMXBs, novae and AGN (especially 
Blazars, which do not show significant intrinsic features). 

6. Conclusion 

We have presented a detailed treatment of the interstellar 
medium towards the nine LMXBs: GS 1826-238, GX 9+9, GX 
339-4, Aql X-l, SAX J1808.4-3658, Ser X-l, 4U 1254-690, 
4U 1636-536, and 4U 1735-444. These bright X-ray sources 
are spread over about 100 degrees in longitude and are located 
near the Galactic plane, which makes them a very good work- 
bench for a study focused on the ISM in the Galactic disk. 

We probed the ISM dust composition, total abundances, and 
abundance gradients. We have shown that along these LOS the 
ISM is composed of a complex mixture of a multi-phase gas, 
dust, and molecules. The analysis of all ionic species of oxygen 
showed that the gas is mostly neutral. The column densities of 



the neutral gas species are in agreement with those found in the 
literature. The column densities that we estimated for the mildly 
ionized gas are generally lower because previous work did not 
take into account that a significant part of the ISM transmission 
may be provided by dust. A significant fraction of oxygen and 
iron is contained in solids, but their exact chemical composi- 
tion is yet to be determined. About 15-25% and 65-90% of the 
total amount of neutral oxygen and iron, respectively, is found 
in dust. This result provides a solution to the problem recently 
raised on the oxygen depletion and total abundance in the inter- 
stellar medium. Interestingly, the dust contribution and the com- 
pounds mixture seem to be consistent along all the lines-of-sight 
(LOS). The ratios between the different ionization states are also 
similar between the LOS. This suggests that on a large scale 
the ISM is chemically homogeneous. Finally, using our X-ray 
data we confirm the abundance gradients in the Galaxy and the 
local proto-Solar abundances as derived at longer wavelengths. 
Through the comparison between our estimates of oxygen abun- 
dance and those obtained by other work on nearby young stars 
and other indicators such as H 11 regions and planetary nebulae, 
we also confirm that the local Galaxy is chemically homoge- 
neous. This work shows that X-ray spectroscopy is a very pow- 
erful method to probe the ISM. 
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Phases 


4U1254-69 


4U 1636-54 


4U 1735^14 


Aql X-l 


GS 1826-24 


GX 339-4 


GX9+9 


SAX J1808.4 


SerX-1 


O i (dust) (n) 


14 


19 


14 


14 


16 


14 


17 


9 


15 


Oi(gas) (a) 


80 


72 


78 


82 


78 


77 


78 


57 


75 


On-v (o) 


5 


7 


6 


3 


5 


7 


5 


17 


9 


Ovi-vm (a) 


1 


2 


2 


1 


1 


2 


<1 


17 


1 


0tot W 


2.10 


2.59 


2.79 


3.31 


2.93 


4.53 


1.66 


1.37 


3.72 



Column ratios are expressed in percentage. ( *' Oxygen total column densities are in linear units of (10 22 m 2 ). 
Table 5. Total abundances of the cold neutral phase (gas + dust). 



Species 4U1254-69 4U 1636-54 4U 1735^44 Aql X-l GS 1826-24 GX 339—4 GX 9+9 SAXJ1808.4 SerX-1 



0/H <n) 


1.12 ± 0.05 


1.16 ± 0.08 


1.14 


±0.09 


1.06 


±0.06 


1.14 


± 0.08 


1.30 


±0.06 


1.27 ±0.09 


1.09 


±0.09 


1.13 ±0.13 


Ne/H (fl) 


0.91 ± 0.04 


1.315+0.07 


1.00 


± 0.07 


0.95 


± 0.03 


1.23 


± 0.06 


1.05 


±0.07 


1.35 ±0.10 


1.37 


±0.08 


1.20 ±0.12 


Mg/H <D) 


1.7 + 0.4 


0.95 ± 0.12 


2.10 


±0.35 


1.00 


± 0.14 


1.47 


± 0.17 


1.39 


± 0.10 


0.8 ± 0.4 


<0.3 


(*) 


1.75 ±0.37 


Fe/H (a) 


0.77 ± 0.36 


1.26 ± 0.21 


0.7 + 


0.3 


1.11 


± 0.12 


1.34 


± 0.17 


1.21 


±0.13 


1.13 ±0.20 


0.77 


±0.23 


1.00 ±0.21 



Abundances ratios are in the linear proto-Solar abundance units oflLodders & Palme (2009). 
For magnesium abundance we obtain only an upper limit. 





Fig. 10. Abundance gradients on a logarithmic scale. The blue dotted lines show the proto-Solar values. 
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